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ABSTRACT 

We extend, for the case of a general scalar potential, the inflaton-graviton 
Feynman rules recently developed by Iliopoulos et al. [[I]]. As an appli- 
cation we compute the leading term, for late co-moving times, of the one 
loop back reaction on the expansion rate for V(<p) = \m 2 ^) 2 . This is ex- 
pressed as the logarithmic time derivative of the scale factor in the coor- 
dinate system for which the expectation value of the metric has the form: 
(0\gnv(t, x)\Q)dx^dx v = —dt 2 + a 2 (t)dx ■ dx. This quantity should be a gauge 
independent observable. Our result for it agrees exactly with that inferred 
from the effect previously computed by Mukhanov et al. |2|, |3|] using canon- 
ical quantization. It is significant that the two calculations were made with 
completely different schemes for fixing the gauge, and that our computa- 
tion was done using the standard formalism of covariant quantization. This 
should settle some of the issues recently raised by Unruh [|J. 
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1 Introduction 



We wish to address a controversy which has arisen in the literature of scalar- 
driven inflation. The dispute concerns the recent claim by Mukhanov, A- 
bramo and Brandenberger 0, [| that infrared modes can generate a significant 
one loop back-reaction which reduces the expansion rate over the course of 
inflation. Unruh M has raised a number of serious questions about their 
methodology and the plausibility of their conclusion. 
We begin by summarizing Unruh's objections: 

1. It is difficult to understand how long wavelength modes can affect the 
local geometry since they should appear spatially constant to a local 
observer. 

2. To leading order in the long wavelength expansion the mode solutions 
are all equivalent, locally, to coordinate transformations which can have 
no effect on local invariants. 

3. The quantization procedure employed by Mukhanov et al. is suspect 
because their dynamical variable is nonzero for only one of the two lead- 
ing long wavelength solutions. Since this dynamical variable possesses 
another, independent solution, the corresponding degree of freedom 
must be unphysical. 

4. What Mukhanov et al. refer to as "gauge independent" quantities are 
really just the local dynamical variables in a particular gauge. 

5. Mukhanov et al. employ an unconventional variation of perturbation 
theory in which the effective stress-energy tensor of the first order equa- 
tions renormalizes the zeroth order stress-energy tensor. 

6. The contributions to the metric at second order — the zero mode of 
which is what Mukhanov et al. computed — depend upon the gauge 
chosen for expressing the first order terms. 

Without wishing to criticize good people who addressed an important 
issue at the extreme limit of their formalism's applicability, we must admit 
to a certain sympathy for Unruh's methodological objections. One of the 
present authors also had difficulty understanding the work of Mukhanov et 
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al. on account of (4-6). [] However, we believe the physics of what they did is 
correct, and that is the point of this paper. 

After fixing notation about the perturbative background in Section 2 we 
comment in Section 3 on the physics of the process and we partially ad- 
dress Unruh's objections (1-3). The remainder of the paper is devoted to 
checking the calculation of Mukhanov et al, in a completely different gauge, 
using the standard formalism of covariant quantum field theory. The Feyn- 
man rules are given in Section 4. These were lifted from a recent paper by 
Iliopoulos, Tomaras, Tsamis and Woodard [jlj, which we have extended so 
that the propagators can be computed (as mode sums) for a general scalar 
potential. Section 5 attaches the external lines (which are retarded propaga- 
tors in Schwinger's formalism |||) needed to convert the amputated 1-point 
functions into the expectation values of the metric and the scalar. We also 
explain how these expectation values are used to compute physical observ- 
ables which measure the cosmological expansion rate and the evolution of 
the scalar. In Section 6 we give the procedure used for isolating the leading 
contribution to each propagator from superadiabatically amplified modes at 
late co-moving times. This is the chief physical approximation of the paper. 
The amputated 1-point functions are computed in Section 7 and processed to 
give the two physical observables. Section 8 summarizes the various results. 

2 The perturbative background 

The system under study is that of general relativity with a general, minimally 
coupled scalar: 

C = j^RV^g - \d^d uVg ^^- g - v(<p)V=g~ ■ (1) 

This section concerns the homogeneous and isotropic backgrounds go and 
ifo about which perturbation theory will be formulated. Three classes of 
identities turn out to be interesting for our purposes: 

1. Those which are exact and valid for any potential V((p); 

1 Other confusing points were the characterization of what is obviously a one loop effect 
as "classical" and the attribution of this effect to an instability in the classical energy 
functional (which is actually stable). 
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2. Those which are valid in the slow roll approximation but still for any 
potential; and 



3. Those which are valid for the slow roll approximation with the potential 



V(<p) = |mV- 



We shall develop them in this order, identifying the point at which each 
further specialization and approximation is made. 

Among the exact identities is the relation between co-moving and confor- 
mal coordinates: 

dsl = -dt 2 + a 2 (t)dx ■ dx = tt 2 (r]) {-drf + dx ■ dx } . (2) 

This implies: 

dt = ttdr] , ao(t) = £l(rj) . (3) 

The Hubble "constant" is the logarithmic co-moving time derivative of the 
background scale factor: 

ao is 

where a dot denotes differentiation with respect to (background) co-moving 
time and a prime stands for differentiation with respect to conformal time. 
Two of Einstein's equations are nontrivial in this background: 

3H 2 = ^ 2 {^o + ^o)} , (5) 
-2H-3H 2 = ^ 2 {^-%)} , (6) 

where k 2 = 167rG is the loop counting parameter of perturbative quantum 
gravity. One can use the two Einstein equations to derive the scalar equation 
of motion: 

V3 + 3Hif + V„(ip Q ) = , (7) 

where V j(p = dV/dip. One can also invert the Einstein equations to solve for 
the Hubble constant and its first derivative: 

H 2 = ^ 2 {^ + ^o)} , (8) 

H = (9) 



Sometimes it is more convenient to write the scalar quantities in terms of 
the Hubble constant and its derivative: 
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0o = —H, (10) 

V( Vo ) = - 2 {H + ZH 2 ). (11) 

At other times one wants to express the scalar quantities using the conformal 
factor ft: 

2 4 1 ( ft" /ft'V. 

= -t? Hf + 2 ill' < 12 » 



ft 



2 i in" /ft'\ 2 ) , x 

^ = ^{n + (n))- (13) 

And the conformal time derivative of the scalar ip' = ft<£> is also useful: 

2 

k 2 \ ft 1 " I ft' 



v'«=-A~ + 4%) ) ■ d4) 



Successful models of inflation require the following two conditions which 
define the slow roll approximation: 

\<p \ < H\tp \ , (15) 
<Pl « ^o). (16) 

It follows that there are two small parameters. Although these are tradition- 
ally expressed as ratios of the potential and its derivatives the more useful 
quantities for our work are ratios of the Hubble constant and its derivatives: 

For models of interest to us the rightmost of these parameters is negligible 
with respect to the leftmost one. We shall also assume that the derivative of 
the scalar is negative: 

(f/ = —QyfH. (18) 

ft 
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The slow roll approximation gives useful expansions for simple calculus 
operations. For example, ratios of derivatives of the field are: 

f = 2tf W ( 1 + J 3 + ...). (20 ) 
Successive partial integration also defines useful slow roll expansions: 

dtH"(l» = jl + (=£U . . . j , (21) 

= 1 1 ™ + ...}. (22) 

J a + lH{a + 2H 2 J V 7 

In discussing the physical significance of their result Mukhanov, Abramo 
and Brandenberger specialized to the simplest potential for chaotic inflation 

§: 

V(<p) -> imV 2 • (23) 

In the slow roll approximation with this potential one can solve explicitly for 
the scalar's evolution: 

2 Tfl 

<fv{t)=ipi--j=-t. (24) 

The interesting geometrical quantities have the following expressions in terms 
of ip (t)- 

H « - l -m\ (25) 
H « -^=«m^ (t). (26) 



Note that H « for this potential, so only one of the slow roll parameters 
is nonzero. The slow roll approximation implies that the initial value of the 
scalar field is much larger than the Planck mass: 

(pi > - • (27) 
Inflation ends in this model when <po(t) ~ hT x . 
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3 Physical comments 



The physical mechanism behind what Mukhanov, Abramo and Branden- 
berger ||, || have found for scalar-driven inflation is roughly the same as 
that studied previously by Tsamis and Woodard in the context of infla- 
tion caused by a bare cosmo logical constant. There is such a simple physical 
model for what is going on that we would be derelict in our duty of explica- 
tion not to present it. Formalists should rest assured that this is merely a 
qualitative description of phenomena whose reality has already been estab- 
lished by computing what should be invariant observables in the standard 
formalism of covariant quantization. 

Owing to the rapid expansion of spacetime and the special properties of 
the dynamical quanta involved^ there is a vast enhancement of the 0-point 
energy which the uncertainty principle requires to be present in each dynam- 
ical degree of freedom. This is the phenomenon of superadiabatic amplifica- 
tion, first studied by Grishchuk |5J. A simple picture for it is that virtual 
pairs with wavelengths comparable to the horizon can become trapped in the 
expansion of spacetime and not be able to recombine. 

Superadiabatic amplification is not a large effect by itself. Although the 
total energy contained in infrared modes increases quite rapidly, the corre- 
sponding expansion of the 3-volume keeps the energy density constant for 
pure de Sitter expansion. f\ For this background it is simple to show that 
there is only about one extra infrared quantum per Hubble volume. The 
interesting, secular effect derives from the gravitational interaction between 
these quanta. As each virtual pair is pulled apart, its gravitational potentials 
fill the intervening space. These remain to add with those of the next pair. 
Even though the 0-point energy stays constant, the induced gravitational 
potential increases. It is the interaction energy between this and the 0-point 
energy, and between the gravitational potentials themselves, which gives the 
effect. 

In the purely gravitational model of Tsamis and Woodard, linearized 

2 These properties are (1) effective massiessness on the Hubble scale and (2) the absence 
of classical conformal invariance. 

3 For certain models of scalar-driven inflation the infrared energy density can grow as 
the scalar rolls. This is what seems to distinguish those scalar potentials for which there 
is a one loop effect from those for which there is not. The two loop effect of pure gravity 
— and presumably also gravity with scalars — does not depend upon such growth. 
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gravitons can only induce gravitational potentials at second order in the weak 
field expansion. Since superadiabatic amplification is a one loop effect this 
means that the secular back-reaction comes at two loop order. When inflation 
is driven by a scalar field its quanta can induce gravitational potentials even 
at linearized order in the weak field expansion. This is why Mukhanov, 
Abramo and Brandenberger were able to follow what is essentially the same 
physical process with a vastly simpler one loop calculation. 

Either way, the effect is to slow inflation because gravity is attractive. 
Since gravity is also a weak interaction, even for GUT scale inflation, the 
process requires an enormous amount of time before it can become significant. 
A direct consequence is that the equation of state of the induced stress tensor 
must be approximately that of negative vacuum energy. To see this consider 
the relation implied by conservation between the induced energy density p(t) 
and the induced pressure pit): 

p = -3H(p + p). (28) 

Since the accumulation of a significant effect requires many Hubble times, 
\p\ <C H\p\, and it must be that pit) nearly cancels pit). 

A sometimes confusing point is that one does not require the complete 
theory of quantum gravity in order to study an infrared process such as this. 
As long as spurious time dependence is not injected through the ultraviolet 
regularization, the late time back-reaction is dominated by ultraviolet finite, 
nonlocal terms whose form is entirely controlled by the low energy limiting 
theory. This theory must be general relativity, with the possible addition 
of some light scalars. It is worth commenting that infrared phenomena can 
always be studied using the low energy effective theory. This is why Bloch 
and Nordsieck |9j were able to resolve the infrared problem of QED before 
the theory's renormalizability was suspected. It is also why Weinberg |T0 



was able to achieve a similar resolution for quantum general relativity with 
zero cosmological constant. And it is why Feinberg and Sucher [1 I] were 
able to compute the long range force due to neutrino exchange using Fermi 



theory. More recently Donoghue [12| has been working along the same lines 
for quantum gravity with zero cosmological constant. 

We emphasize that the process is causal, in spite of its close association 
with modes whose wavelengths have redshifted beyond the horizon. This 
emerges most clearly in the two loop computation of Tsamis and Woodard 
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where the effect derives from integrating interaction vertices over the past 
lightcone of the point at which the expansion rate is being measured. Be- 
cause gravitons are massless these interactions superpose coherently. Because 
gravitons are not conformally invariant they reflect the enormous physical 
volume of the past lightcone rather than its minuscule conformal volume. 
The growth in the back-reaction is directly attributable to the fact that the 
invariant volume of the past lightcone increases without bound as one ob- 
serves at later and later times. 

Causality is also built into the work of Mukhanov, Abramo and Branden- 
berger through their use of the Heisenberg field operators. The equations of 
motion for these are simply operator realizations of the causal field equations 
of classical general relativity. In a local gauge one can express an operator 
at the spacetime point (t, x) entirely in terms of the operators and their time 
derivatives on that part of the initial value surface which lies on or within the 
past lightcone of (t, x). There might be some dispute about what this means 
nonperturbatively, where the quantum metric can have a significant impact 
on the lightcone, but it makes perfect sense in the perturbative regime under 
study. The time dependence Mukhanov, Abramo and Brandenberger obtain 
derives in part from the continual redshift of new modes from ultraviolet to 
infrared but mostly from the growth of the infrared mode functions which 
can occur in some (but not all) inflationary backgrounds. 

We turn now to issues (1-3) listed in our Introduction. Regarding the 
first objection, it is relevant to note that while long wavelength modes indeed 
appear spatially constant to a local observer, so too does the cosmological 
expansion rate. Therefore the causative agent and its purported effect are 
commensurate. 

Viewed from the perspective of obtaining a long period of inflation it is 
rather local, short wavelength phenomena that ought to be regarded with 
suspicion. Without severe fine tuning the natural duration of any process 
mediated by short wavelength quanta must be the Hubble time or less. It 
seems reasonable to conclude that a mechanism for screening the cosmological 
constant must also end inflation if nothing else does the job first. But inflation 
has to persist for many Hubble times in order to explain the large scale 
smoothness of the observed universe. Note that infrared, long wavelength 
phenomena can require much longer to produce a significant effect because 
they can act by coherently superposing an inherently weak interaction over 
the past lightcone. It is only an enormous expansion of the invariant volume 
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contained in the past lightcone that can compensate for the weakness of 
gravitational self-interactions. 

Finally, one must distinguish between a local observation of the cosmo- 
logical expansion rate, and the local expansion rate that would be produced 
by spatial inhomogeneities in the vacuum energy. Many people believe that 
whatever is suppressing the former must also suppress the latter. We do not 
share this view. Experiment is sadly unable to decide the matter but it seems 
to us that a local fluctuation which created a large enough region of negative 
p + 3p should result in that region beginning to undergo inflation. We believe 
that known physical principles already suffice to explain why such fluctua- 
tions are rare in the observed universe |13|. Were it otherwise one would not 
be able to make conventional models of inflation agree with observation by 
the unaesthetic device of fine tuning the bare cosmological constant which 
is, be it noted, spatially homogeneous. 

The relevant point about the second objection is that the "stress tensor" 
of gravitational perturbation theory is not an invariant, or even a scalar. So 
the fact that infrared mode solutions look, to leading order for small wave 
number, like coordinate transformations does not mean they necessarily have 
no effect. Superadiabatic amplification allows these modes to carry nonzero 
energy and pressure in spite of their extreme redshift. The proper way to 
determine their effect is by computing the metric 's response at second order 
and then addressing the gauge issue of what this response means physically. 
We believe this is what Mukhanov, Abramo and Brandenberger did, although 
perhaps not as transparently as one might wish. To check their result we 
made what ought to be the same computation in a completely different gauge 
and using the standard formalism of covariant quantization, and we got the 
same answer. 

This is the right point to comment on the gauge issue, which was also 
raised extensively by Unruh. What both we and Mukhanov et al. computed 
was the expectation value of the metric in the presence of a particular state 
and in a fixed gauge. There is no doubt that this quantity depends upon 
the gauge in which the computation was done. It is important to realize 
that a quantity is not automatically devoid of physical import by virtue of 
being gauge dependent. It can still contain useful physical information which 
can be separated from the unphysical, gauge dependent part. Examples of 
this abound in quantum field theory. The most straightforward is the way 
in which gauge dependent Green's functions can be processed, using LSZ 
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reduction, to give gauge independent, on-shell scattering amplitudes. (This 
is discussed in any standard text on quantum field theory, for example that 
by Peskin and Schroeder JT3).) One does not even have to consider products 
of field operators. There is an elegant formalism, due to DeWitt [|15[], in which 



the S-matrix is obtained from the in-out matrix element of the dynamical 
variable in the presence of a general scattering state. 

So the expectation value of the metric contains valid physical information; 
the question is how to extract it. Our technique exploits the special property 
of the initial state of being homogeneous and isotropic. This means that a 
co-moving coordinate system exists for which: 

(0\g lw (t,x)dx IM dx v \Q) = -dt + a 2 (t)dx-dx . (29) 

Our observable is the logarithmic time derivative of the scale factor in this 
coordinate system: 

One can investigate how this quantity changes under a variation of the gauge 
fixing functional and the result is that it does not change 0. This would 
seem to be the analog of De Witt's theorem about the gauge independence of 
the on-shell S-matrix. Of course the absence of gauge dependence does not 
automatically endow a quantity with physical import. We interpret H e{i (t) 
as the expansion rate a local observer would measure in the presence of state 
|0). It certainly has this meaning in the classical limit but we are willing to 
entertain dissident views. 

We come finally to Unruh's doubts about the formalism of Mukhanov, 



Feldman and Brandenberger ||16|| . His argument is based on the long wave- 
length solution [|J he found for the linearized Newtonian potential, which we 
shall call n(x). In our notation this quantity corresponds to the following 
invariant element: 

g^dx^dx" = Q 2 {-(1 + 2Q- 1 n)dr] 2 + (1 - 2Q- 1 n)dx ■ dx\ . (31) 

When the linearized goi equations are used to eliminate the scalar field the 
linearized g 00 equation becomes: 

n " _ V 2 n - 2^n' + — - 2— n = . 32 
tp \\l ir J 
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In the limit that the V 2 term can be neglected Unruh obtains the following 
independent solutions: 

"i = ft-^/^Ajtffl)), (33) 
- = (34) 

Unruh's problem concerns what happens when n\ and n 2 are substituted 
into the dynamical variable used by Mukhanov, Feldman and Brandenberger: 

4 i f , / n" n'\ } 

Although n\ produces a reasonable function: 

4 1 f QQ" A ft 2 

the n 2 solution gives v = 0! This is disturbing because the formalism of 
Mukhanov, Feldman and Brandenberger quantizes v as a scalar field which 
obeys the following second order equation: 

v » _ V \ _ < v = o . (37) 

In the limit that the V 2 term can be neglected one finds that v — v\ is indeed 
a solution, as is: 

Since this second solution does not correspond to any combination of Unruh's 
two long wavelength solutions he concludes that it must be unphysical and 
that the formalism is therefore suspect. 

In fact neither the v-i solution nor the n 2 solution is unphysical, they 
simply correspond to different orders in the long wavelength expansion. At 
fixed, nonzero wave number k one can express the two independent solutions 
for the Newtonian potential as power series in k 2 . The zeroth terms in these 
two series are Unruh's solutions, 711(77) and 712(77). However, one should really 
include some higher order terms as well since the physical relevance of the 
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solutions is for modes with small but nonzero wave number.^ The first order 
correction to n 2 can be expressed using an advanced Green's function: 

4A; 2 f® ^2(^7) 
N 2 (rj,k) = n 2 (v) + — / drj {wi (77)712(77) - n 2 (v)m(v)} , 77=^2 +0(k) . 

(39) 

When it is substituted for n(x) in ( p5|) the result is: 



4fr 2 

^U=v 2 = — v 2 ( V ) + 0(k 4 ). (40) 

So the solution sets of the two variables are in one-to-one correspondence 
and there is no obvious problem with the formalism of Mukhanov, Feldman 
and Brandenberger. 

At the price of specializing to power law inflation one can even see how 
the n and v solution sets relate to all orders. For Q = {rjo/rjY'^ 1 ' the 
two mode solutions for the Newtonian potential are proportional to Bessel 
functions of order /x=| + l/(s — 1). Unruh's solutions are the zeroth order 
terms in the power series expansions of the following: 

N^k) = (41) 
N 2 (v,k) = 8 —r(l + l i)(!^) "-^Uk V ). (42) 



/VoV 

One can easily verify that the solutions to ([57]) which make contact with v\ 
and v 2 are: 

V lM = llt^l^f 1 (43) 

K y/S \ 2 J \/ 7/0 

V 2 (v,k) = -^V*T(1+A*)fe) " ■ (44) 



4 There are a lot more of these. The number of zero modes is constant in time whereas 
the inflationary redshift eventually makes the physical wave number (f7 _1 fc) of any mode 
small. Further, as Unruh pointed out, the k = system is degenerate in that his two 
solutions become unphysical on account of being exactly coordinate transformations. The 
one physical solution is paradoxically absent from the k 7^ system. It appears because 
the goi equation is automatically satisfied for k = and therefore fails to relate the zero 
modes of the Newtonian potential and the scalar 01 . 
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For power law inflation the relation fl3"5| ) between n and v is recognizable 
as the recursion relation which produces — J^+i from J M (and +J_ /Lt _i from 

v^~^y/v\[y/m)'-^[yM] > («) 
and the V^'s descend from the AVs as follows: 

Vjfo,*) = -l^VvllVvN^W-^m^k)}) , (46) 
^fo.*) = ^x—^V^I^^,*)]'-^^,*)]} -(47) 

4 Feynman rules 

The purpose of this section is to give the Feynman rules for the general 
inflaton-graviton action ([!]). We have mostly borrowed these from a recent 
paper by Iliopoulos, Tomaras, Tsamis and Woodard H. The one exception 
concerns the issue of mixing between the scalar and the 00 component of the 
graviton field. Iliopoulos et al. were only able to solve the system for a class 
of backgrounds including those of power law inflation, but not the power law 
backgrounds typical of chaotic inflation. We have achieved a general solution. 
One should also note that all formulae given in this section are exact. We 
have made neither the slow roll approximation nor have we specialized to the 
case of a quadratic potential. Of course that will be necessary in order to 
convert the formal mode sums into explicit results, but the task of making 
these approximations has been postponed to the end of Section 5. 

Our quantum fields are the scalar and the conformally rescaled pseudo- 
graviton 

<p = (f + <f), (48) 
= Q 2 (tj^ + m/v) = Q 2 g^ . (49) 

It should be noted that cosmologists typically restrict the word "graviton" to 
that part of the metric which interpolates dynamical spin two quanta at lin- 
earized order. Adhering to this convention would be terrifically cumbersome 
in the context of BRS quantization beyond linearized order. Our "pseudo- 
graviton" also includes degrees of freedom which are constrained or pure 
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gauge. This is the standard usage in particle theory, cf. the "photon" field, 
propagator and interactions of QED and the "gluon" field, propagator and 
interactions of QCD We shall try to avoid misunderstandings, without 



over-burdening the notation, by following the convention of Iliopoulos et al. 
who parenthesized the word "general" before "pseudo-graviton." 

As usual, (general) pseudo-graviton indices are raised and lowered with 
the Lorentz metric. After many tedious partial integrations the invariant 
Lagrangian can be written as a total derivative plus the following: 



- \\F^9 pa r v ^,^ v a {^), a - nVo^y^ 

where a comma denotes differentiation. Gauge fixing is accomplished by 
adding a gauge fixing functional and the corresponding ghost action to obtain 
the BRS Lagrangian: 

£brs = Anv - \rTFvFu - n W^SFp . (51) 

The symbol SF^ represents the variation of the gauge fixing functional under 
an infinitesimal diffeomorphism parameterized by the ghost field u^. We will 
follow Iliopoulos et al. in our choice of gauge fixing functional: 



1 Q' 

F» = ft ( - 2^ - 2 q-^o + ^o^o0 ) . (52) 



A great advantage of this gauge is that it decouples the tensor structure of the 
propagators from their dependence on spacetime. The propagator becomes 
a small number of constant tensors multiplying only three different types of 
mode sums. Another advantage is that the limit Q — > 1 takes this gauge 
to one of the standard gauges of flat space, which often provides a useful 
correspondence check. 
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With a few more partial integrations the terms quadratic in the various 
quantum fields can be reduced to the following form: 



1 ,rl" 

ft* + i K 2 ,2 <£g_ 



^ D A - tJ v D B 



(53) 



A number of pieces of notation require explanation. The differential operators 
Da and D# are: 



D A = n 
D B = n 



d 2 + 
d 2 + 



in- 



n- 



i 



n 



(54) 
(55) 



where d 2 = rj^d^ is the d'Alembertian in conformal coordinates. It is 
worth commenting that Da is the kinetic operator for a massless, minimally 
coupled scalar. The kinetic operator for the (general) pseudo-graviton is: 



D pa 



Da - t^5 u) (p t^D B + t^tn^B . 

(56) 



2 n » 4 «#~ > 2 
Parenthesized indices are symmetrized, the symbol denotes: 



(57) 



and a bar above a Lorentz metric or a Kronecker delta symbol means that 
the zero component is projected out: 



7/ 



5; + . 



(58) 



The quadratic Lagrangian involves two sorts of mixing: that between the 
spatial trace and ^oo and that between ^oo and the scalar 0. The first can 
be removed by the following simple field redefinition: 



Ipij = Cij + %CoO 



IpOi = Cl 



" 



^00 = Cl 



00 



(59) 
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where small Latin letters denote spatial indices. In these variables the 
quadratic part of the Lagrangian becomes: 



(2) 
BRS 



I in'" 
+ £ K 2 /2 fQ_ 



(60) 



For a general Oo(t) there is no local change of variables which removes the 
mixing between (oo and <p off shell. However, it is easy to diagonalize the lin- 
earized field equations which determine the on shell mode solutions. Canon- 
ical quantization of the linearized theory can then be invoked to expand the 
original quantum fields in terms of creation and annihilation operators. It is 
straightforward to use these expansions to express the propagators as mode 
sums. At the cost of some mixed propagators we will eventually obtain a 
complete expression of the Feynman rules in terms of the field variables ipfiu, 
ujy and 0. 

It is simplest to absorb a factor of Q into Coo an d <$'■ 



z = fXoo , f = ^4> 
The linearized equations for the z—f system are: 



d 2 + £«V 2 ) z ~ K( Pof 



-Ktf%z+ id 2 + -«Vd 2 



f 








(61) 

(62) 
(63) 



Differentiating (p^) and adding it to ^Kip' Q times 
diagonalized field equations: 



gives the first of our 



2,./2 



d 2 + 7 *Vo 



. 



(64) 



The second diagonalized field equation comes from differentiating (|63|) and 
adding it to \^' times ( |62| ) minus y^'/Vo times (|63|) : 



V + i.V 2 -^ + 2^' 



(65) 
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The preceding discussion implies that the diagonal variables are: 



x(rj, x) = z'(r], x) + -Kcp' (r])f(r], x) , 
y(r),x) = ^Kip' (r])z(r],x) + f'(r],x) - 



(66) 
(67) 



Since conformal time derivatives appear in the transformation its inverse 
cannot be local in time for the off shell fields. However, by using the linearized 
field equations one can obtain the following expressions for the conformal 
time derivatives of x and y: 



x 



y 



(V 2 + ±«V 2 ) Z + \K^f - l -K^f , 



(68) 



1 «JJ - Wo* - 7Tf' + ( V2 + Wo 2 + %) f - (69) 



2 ru 2 " " y?b 
Eliminating z' and /' gives the following on shell inverse transformation: 

, 1 



z 



1 

V2 



x - 2 K ^y 



f = 



i 



--K<p' x + y' + ^-y 



(70) 
(71) 



We stress that since the linearized field equations have been used these rela- 
tions apply only to the on shell mode solutions, not to the off shell fields. 

The mode equations for all the fields — including flQj, fl(oi, fiu7 M and 
VLuj u — can be given a simple, unified treatment. There are three types of 
modes which we shall call A, B and C. They are defined as the plane wave 
solutions annihilated by the following differential operator: 



Vj EE d 2 + 



0i 



(72) 



where k = 

e A = n 



and the various 0i(r)) 's are: 
, 9 B ee Q 1 , 9 C 



2 n' 



H 



k VL 2 ip' 



= a r 



(73) 
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Briefly, the spatial polarizations — 0^-, QuJ 1 and Vtujj — are comprised of 
A modes, the mixed polarizations — ^Coi, and Quq — are made of 
B modes, as is the diagonal variable x, and the other diagonal variable y 
consists of C modes. Because quantization was accomplished by adding a 
gauge fixing term most of the linearized fields harbor unphysical quanta. 
Physical gravitons are A modes that reside in Qj] the physical scalar is a C 
mode in y. 

We will return in Section 6 to the problem of obtaining useful approx- 
imations for the mode functions but we proceed, for now, as though they 
are known. We define Qi(rj,k) as descending by perturbative iteration (ex- 
plained in Section 6) from the pure negative frequency solution for wave 
number k = \\k\\. We also assume it has been canonically normalized: 

Qi(v, k)Q*/(v, k) - Q' T {v, k)Q* T (v, k)=i. (74) 

From canonically quantizing the quadratic action (|60|) one finds that the 
fields z(r), x) and z'(r), x) form a conjugate pair. The same is true for f(jj, x) 
and f'(f],x), so the only nonzero equal-time commutators involving these 
fields are: 

[z(rj, x), z'(r], y)\ = id 3 (x - y ) = [f(r], x), f'(r], y)\ . (75) 



From their definitions (|66| - p7T) and the on shell relations (p8|- p9|) one can easily 
check that the only nonzero equal-time commutators in the x-y sector are: 

[x(?7, x),x'(r], y)\ = -iV 2 5 3 (x - y) = [y(rj, x),y'(rj, y)\ . (76) 

We can realize these commutation relations with conventionally normalized 
creation and annihilation operators (X, JT' and (Y,Y^): 

~X(k),X\p)} = {2-nf5\k-p) = [Y(k),Y\p)} . (77) 

Since T>Bx(r],x) = we expand x using B modes: 

x(v,x) = J 0- 3 ke rk - £ {x(k)Q B ( v ,k) +X\k)Q* B (v, k)} ■ (78) 
Since Vcy{fj,x) = we expand y using C modes: 

V(V, x) = J ^0r 3 ke^ { Y (k)Qc(v, k) + Y\k)Q* c (v, k)} . (79) 
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The on shell transformations ( |7D| - [7TD allow us finally to give operator expan- 
sions for z and /: 



z(r],x) 



ik L 



d 3 k e 



(2vr) 3 k 



{-X(k)Q' B ( V ,k) 



d 3 k 



ik-x 



(2tt) 3 k 



+Y(k)-K(p' (r])Qc(v,k) +c.c. 
X(k)^' ( v )Q B ( v ,k) 



(80) 



-Y(k) 



d?7 ^,(77) 



Qc(?7, k) + c.c. 



(81) 



The rest is a standard exercise in free field theory. We choose the state 
|0) to obey: 

X(£)|0) = = Y(£)|0) . (82) 

The various propagators can be most conveniently expressed in terms of the 
mode sum iSi(x; x'): 5 



iSj(x;x') = - — 



d 3 k , 



ik-Ax — ek 



(2tt)* 



{OiArfiQfakWjWk) 



2tt 2 Jo 



dk 



+9(-Ar ] )Q*( V ,k)Q I ( V , ,k)} , (83) 
sm(kAx) ^_ ek 



kAx 



-e-^ie^Qj^Q^'.k) 



+e(-A V )Q*(v,k)Qi(v',k)} 
where we define the following conformal coordinate differences: 
Ar) = rj — 7]' , Ax = x — x 1 , Ax 



34) 



(85) 



Expanding the various time-ordered products and exploiting fl32| ) leads to 
the following expressions for the propagators: 

(0\T{z( V ,x)z( V ',x')}\0) = ^-^5 B (x;x') 



3 Note that we have introduced a time independent, ultraviolet convergence factor of 
This corresponds to an exponential mode cutoff on the initial value surface. 



-eft 
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+ -K 2 <p' (r))<p' (r)')i5 c (x;x') , (86) 



(0\T{z(r),x)f(r)',x')}\0) 



d_l 

"drj2 



K<p' (r)')i5 B (x;x') 



drf (f' 



{0\T{f(r,,g)f(rf,^)}\0) = -«VoMvW)*M^') 



+ 



dr] iff Q 



dr}' ip! Q 



i5 c (x; x 1 ) , (87) 



i5 c (x;x') . (88) 



Now we restore the factors of Q, 1 . These can be used to convert the 
conformal time derivatives to co-moving time: 



The three z-f propagators become: 



Q (p' = ipo ■ 



(89) 



iA a (x;x') = (0|T{^oo(»7,^oo(»/,^)}|0) 
(9 (9 1 

iA^x-x') = (O|T{^oo(^x)0(V,f')}|O) 



9 1 1 

= —^^0o(t')iS B (x; x') - -mpo(t) 

tA,(x;x') = (O|T{0(77,x)0(7/,f)}|O) 
= -K 2 (po(t)(p (t')iS B (x;x') 



of (f 



(90) 

i^fox 7 ) , (91) 



+ 



9 +H+^ 



dt 



iSc(x-x') . (92) 



The unmixed propagators can also be represented as mode sums: 

iA I (x;x) ee fr^frVX-v 2 )^/^) • 

The (general) pseudo-graviton propagator is: 



(93) 



(0|T{Vv(x)V>(<>} |0> = tA A (x;x')2 



Vfi(pVa)u ~ VfiuV 



pa 



4iA B (x; x')t^rj u u p t a) + iA a (:r; a/) rj + r? p(T + t p t a . (94) 
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The other propagators are: 

(O|T{Vv(z)0(z')}|O> = i^{x-x' 

(0\T{<p(x)<P(x')}\0) = iA^x'), 

(0\T{u,(x)uJ v (x')}\0) = iA A (x;x')rJ^ 



V + V"] ' ( 95 ) 

(96) 

iA B (x;x')tpt u . (97) 



All interactions between (general) pseudo-gravitons and scalars can be ob- 
tained by expanding the invariant Lagrangian Q5DD using the following iden- 
tities: 



1 + + KM -V> - + • • • • 



(98) 
(99) 



Interactions involving co> M and u v can be read from the ghost Lagrangian, 
which we have simplified by neglecting some total derivatives: 



C 



(100) 



ZJ P 



5^ D A - t^B B 



I 2i 



91 



(101) 



5 Attaching external lines 

The purpose of this section is to explain how we pass from the amputated 
1-point functions which are actually computed to physical observables. We 
begin by expressing the effective Hubble constant in terms of the (nonampu- 
tated) 1-point function. The rest of the section is devoted to the procedure 
for attaching the retarded propagators, needed in the Schwinger formalism 
|J, to convert amputated 1-point functions into their nonamputated cog- 
nates. Although an exact solution is obtained we specialize it, at the very 
end of the section, to leading order in the slow roll approximation. 

Both the initial state and the evolution equations are homogeneous and 
isotropic. It follows that the expectation values of the (general) pseudo- 
graviton field and the scalar can be expressed in terms of three functions of 
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rj. It will simplify some of the later formulae if we choose to think of these 
quantities as functions of the co-moving time t of the perturbative back- 
ground, even though we still are expressing them in conformal coordinates: 6 

(0|«WM|0) = A(t)fJ^ + C(t)t^U, (102) 
(0\K(f)(r),x)\0) = D{t) . (103) 

Note that we cannot assume C = —A since the expectation value may not be 
conformal in the perturbative coordinate system. None of these quantities 
is itself physical but they can be combined to produce observables. We first 
construct the invariant element to infer the true scale factor and co-moving 
time 1 of the expectation value of the metric: 

- dt + a 2 (t)dx -dx = fi 2 {- (1 - C) drf + (1 + A) dx ■ dx} , (104) 

= -[1 - C(t)]dt 2 + a 2 (t)[l + A(t)]dx ■ dx . (105) 

One physical observable is the effective Hubble constant expressed as a func- 
tion of the co-moving time t: 

H M t) = |l,KDl = [H(t) + i^L} . (106) 

If the scalar can be measured then its expectation value is also an observable 
when expressed as a function of the co-moving time t. We shall call this 
variable $(t): 

$(t) = <p (t) + -D(t) . (107) 

K 

What we actually compute are not the expectation values of Kip^ and K(f) 
but rather the amputated expectation values with the external propagators 
removed. We will use Greek letters to denote the three functions of t which 
describe these amputated quantities: 

«g)y + 7(0^ = D^ PC J (0 1 Kipp* M|0) 

6 Recall that the relation between t and r\ is: 

dt = fl(r])dr] <=^> = drj 

a (t) 
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-kQVoV"(°M^)I > 
S(t) = -^VoW(0|^ pCT (r7,f)|0) 



108) 



V 4 <PoJ 



Contracting with the kinetic operator 
ponents gives three relations: 



a 

7 
8 



—V A (A-C), 

^B A (A -C) + D B C - KfiVoZ> , 



d 2 



drj 2 



fi(O|K0(i7,x)|O) . (109) 
and isolating distinct tensor com- 

(110) 

(111) 

VLD . (112) 



1 in'" 

4" ^° + ^ 



Since it is from A and C that H ef[ is constructed, we must invert these 
relations. 

We employ the Schwinger formalism || in order to get true expectation 
values rather than in-out matrix elements. An important feature of this 
formalism is that external legs are retarded propagators. This means that the 
coupled differential equations in (|110H112|) must be inverted using retarded 
boundary conditions: 



= A(0) = C(0) = D(Q) , 
= A(0) = (7(0) = D{Q) , 



(113) 
(114) 



Now it happens that every differential equation we have to solve can be cast 
in the form: 

' + ?1 /(*) = 9(t) ■ 



drf 







(115) 



This is fortunate because the retarded solution can be simply expressed as a 
double integral: 

f(t) =V-\g) = -9( V ) f dt^it^O- 2 ^) P dha^ihWriMh) . 



(116) 

It simplifies the algebra somewhat to multiply the amputated quantities 
by and their unamputated descendants by Q. We denote the rescaled 
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variables by a tilde: 

A = QA , C = QC 

5 = £l~ a , 7 = f2 _1 7 , 

In this notation the equations we must invert are: 



D = QD , 



(117) 
(118) 



a 
1 
8 



--V A {A-C), 
-V a {A-C)+V b C-k^D 



d 2 



drf 4 ip 



D 



(119) 
(120) 
(121) 



where V>a and T>b have the form ( |115|) with 9 a = ^ and 6b = Q 1 - Equation 
( [TT9D implies: 

A = C + V2 l (-45) . (122) 
Substituting this into ( |120| ) gives: 



35 + 7 = V B C - Kip'^D 



(123) 



which, with ( |121| ), is similar to the coupled z—f system of Section 4. Paral- 
leling the analysis of that section we differentiate ( [L23| ) and add it to 
times (|121| ) to obtain: 



V u ( C + = 35' + f + . 



(124) 



Differentiating (|121|) and adding it to ^nf' Q times ( |123|) minus <p'o/<p' Q times 
( PID gives: 



V C (\k^C + D>- &d) = i (35 + 7) + ? - 4^ , (125) 

where £> c has the form ( [H5|) with C = -2Kr x Qr 2 Q! f iff^ 

Of course we can invert the differential operators in the last two equations: 



a + 1 -kv' q d = v B l (35' + 7' + l^' s ) 



(126) 



+ = P c 1 fW (35 + 7)+5'-4^ , (127) 
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but this still leaves derivatives on C and D. These derivatives cannot be 
removed as we did in Section 4 because the Laplacian vanishes for spatially 
homogeneous functions. What we must do instead is to divide ( |123| ) by ip' 
and add it to 2<// '/<^q 2 times (|126|) : 

{ 4) = ^ + (35' + ~i + . (128) 



Dividing (|121|) by y?o and adding it to 2<p'q/<p' times ( [L27|) gives a similar 
relation for D: 

f = £ + 2-^Z# f \ Vo(35 + 7) + ? - 4*) • (129) 

Putting everything together gives the following solution for the unampu- 
tated coefficient functions: 

c = ^,f^ +2 <£^L, +f+ i ^\ (131) 



11 { if (p' V 2 

D = #^ 1 {4 + 24^ 1 f^(35 + 7) + ?-^)| . (132) 

Recall that a = 7 = f2 _1 7 and 5 = The various inverse 

differential operators are defined by (|116|) with the following assignments for 
9(rj): 

9 A = n , e B = Q- 1 , e c = --^rr. (133) 

For the potential V((p) = \m 2 (p 2 it happens that the leading order results 
for the amputated 1-point functions consist of sums of terms with the general 
form: 

a{t) = a N H N (t)4(t) + . . . , (134) 
7 (t) = lN H N (t)at(t) + . . • , (135) 



Sit) = H X U^^-\ .... (136) 
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The coefficients oln, 7n and 5n are constants. When the slow roll expansions 
of Section 2 are applied to the various integrations and differentiations in our 
formulae for A, C and D, the following leading order results emerge: 



A N (t) = 
C N (t) = 

D N (t) = 



+ ... 



3N { H\t) j \-H / 
3a N + ^ N \ (3a N + 7^ + 5 N " 



'3a N + 7at + 5 



N 



3N 
'iff - 



fff - H (t) 



H (t) 



H*(t) 
H 



3N 



H\t) 



+ . . 



(137) 
+ • (138) 
. (139) 



Here Hj = ff (0) is the Hubble constant at the beginning of inflation. These 
results imply the following leading order shift in co-moving time: 



(t-t) 



N 



3 ( 3oln + 7jv + 8n 
2 aN -[ 3N 



'iff 



N Hf 



N-l 



N — 1 H 



+ 



\ R2 
H N ~ 2 \ I H 



N-l 



+ ... 



-H, 



To leading order the proportional shift in the two observables is: 



ffeff ~ H \ 

H J 



N 



-a N H + 



3a N + 7at 



N 



'3a N + -f N + 5. 



N 



6(N-1) 



H 



N-l 



.N-l 

R N-l\ 



H 



+ ... , 



(140) 



(141) 



N 



3a N + 7at 



N H 



N-l 



iff 



ff 2 N-l H 



+ 



H 



N-2 



N-l 



'3a N + j N + 5 



N 



6(N-1) 



H -l H N-l _ jj 



N-2 



(142) 



To obtain the full shift one sums the contributions for various different values 
of N. 
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6 Infrared parts of propagators 



This section deals with a very important omission in Section 4. Although we 
were able there to express the various propagators as mode sums for an ar- 
bitrary scalar potential, we do not possess the corresponding mode functions 
Qi(j],k) for a general potential. This is a standard problem in the theory 
of cosmological perturbations WB, [L7j and we solve it in the standard way: 



by developing series solutions for the ultraviolet (early time) and infrared 
(late time) regimes. The normalization for the ultraviolet expansion derives 
from the flat space limit. We normalize the infrared expansion by matching 
its leading term with that of the ultraviolet expansion at the time when the 
physical wavelength of each mode is just redshifting beyond the Hubble ra- 
dius. (This is the chief approximation of the paper.) One then defines the 
"infrared part" of each propagator as that obtained from the leading order 
term of the infrared expansion. We report explicit results to leading order in 
the slow roll approximation. 

Recall from Section 4 that we have three kinds of plane wave mode solu- 
tions Qi(fj, k). They obey the equation: 

ViQiiv, k) = - { ^ + - Qiiv, k) = 0, (143) 

where the ^/(//)'s are: 

e A = n , e B = n- x , e c = --^ T = a l ^=. (144) 

« " Vo J-H 



Of course (|143|) does not completely define the modes because there are two 
linearly independent solutions. We define Qi(r),k) as the solution of ( |1 13| ) 
which is canonically normalized: 

Q( V , k)Q*'(n, k) - Q'(r], k)Q*(r], k) —i , (145) 

and descends by perturbative iteration from the negative frequency solution 
of the far ultraviolet. 

The far ultraviolet is defined by k 2 ^> 6" /6. At fixed k this condition will 
also be realized, in all models of inflation, as the conformal time approaches 
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negative infinity. In the ultraviolet regime we build up normalized solutions 
by iterating the following equation: 

Qfa, k) = ^e- ikv + /£ d*l\ sin [k(rj - 77)] jMqM k) . (146) 

The result is a series in inverse powers of k. These solutions are obviously 
negative frequency in the far ultraviolet. Their Wronskian ( [L45|) is constant 
as a simple consequence of the mode equation ( 143 ) while its actual value 
derives from the fact that 9"/9 vanishes as the conformal time approaches 
negative infinity. 

The far infrared is defined by k 2 <C 9"/9. At fixed k this condition will 
also be realized, in all models of inflation, as the conformal time approaches 
zero from below. One can find explicit solutions in the limit that the k 2 term 
is neglected. The first one has the same form for I = A,B, C: 

Q 10tI {r]) = 9j( V ) . (147) 

The second is an integral whose convergence (for models of inflation) requires 
different limits for I = A: 

ro dr} 



QioAv) = -0i(v) I ^z=y > ( 148 ) 

and I = B,C: 

Qw,i(v) = di(v) £ i (I = B,C). (149) 

When k 2 is small but not zero one can build up solutions which descend from 
Qio,i{v) by iterating with the appropriate Green's function: 

Qij(r],k) = Q i0 ,i(ri) ~k 2 df] G app (r],rj)Q iJ (r],k) (i = 1,2) , (150) 

J — 00 

This obviously gives a series of increasing powers of k 2 . Here the "appro- 
priate" Green's function is chosen to make the integral converge. The four 
possibilities are: 

G adv (r], rf) = +8(7) ~ v)Qw(v)Q2o(ji) ~ ~ v)Q2o(v)Qw(jj) , (151) 
£12(77,77) = -9(V ~v)Qio(v)Q2o(v)-9(t]- ^20(^0(7]) , (152) 
^21(77,77) = +^(77 - 77)g 2 o(77)g 1 o(77) + ^(77 - 77)g 10 (77)g 2 o(77) , (153) 
G re t(77, 77) = +9(v-v)Q2o(v)Qio(j])-0(v-v)QMQ2o(j]) , (154) 
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and it should be noted that one may have to switch from one to another 
midway through the iteration process. 

Since the full infrared solutions, Qij(rj, k) and Q2,i(Vi k), span the space 
of solutions to V I = 0, it must be possible to express the ultraviolet solutions 
as linear combinations: 

Qi(r), k) = qiQ 1: i(r], k) + q 2 Q2,i(Vi k ) ■ (155) 

If we had the full solutions it would be straightforward to determine the 
combination coefficients: 

Qz/jh k )Qi(v, k ) - Q2,i(v, k)Q'j(r), k) 

-Q'i,i(v, k]QAjh k) + Qyfa, fc)Q' f fa k) 
q2 " Q l Mk)Q , 2Av,k)-Q2,i{v,k)Q'iMk)' [ ' 

where any conformal time rj could be chosen. 

For most backgrounds we do not possess the full solutions — either in 
the ultraviolet or the infrared. However, it happens that one of the zeroth 
order infrared solutions — either Qwjiv) or Q2o,i(v) — dominates the other 
and all corrections as the conformal time approaches zero from below. The 



standard approximation [16, 17 1 is to match this solution with the zeroth 
order ultraviolet solution at the horizon crossing time 77*, whose defining 
condition is: 

k = H*n(rj*) . (158) 

Then the behavior of the modes in the far infrared can be approximated as 
follows: 

Q.M^^f^, (IM) 
QioAV*) V2k 



where i is either 1 or 2, depending upon which of the zeroth order infrared 
solutions dominates for the I mode. 

For J = A it is Qio,a(v) = ^iv) that dominates at late times. We can 
therefore write: 

tt -ikij* 

q a ^ k ) — Q(r,) ■ -± (160) 
k V2k 
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For I = B Q2o,b = 1 becomes irrelevant at late times. The dominant 
solution is: 

Q20AV) = A" 1 f ^ dvQ 2 (rj) = -L f ^ dt'a (t') « -L . (161) 
We can therefore approximate the B modes as follows: 

Q ^- k ^WY H '^- (162) 



The case of I = C requires a more extensive analysis. The first solution is 
down by an inverse scale factor but enhanced by the inverse of a slow roll 
parameter: 

Q\o,c = a o 1 / — = ■ (163) 
-H 



The second solution has to be re-expressed several times before it can be 
recognized as the dominant one: 

Q 2 „. c = a^4=J[jt% >{ t^y^ , (164) 



1 d 



jj dt (^Clg(t) J— 00 



[ dt'a (t')\ , (165) 

J — 00 



-H 



H 2 



(166) 



After any significant amount of inflation the inverse scale factor is much 
smaller than —H/H 2 , so we can approximate the C modes as: 



Qciv, k) — ► ±— r=~/ W ■ ( 16T ) 



Although the approximations ( p.60|J162| , |167| ) we have just made may seem 
grotesque they are intimately related to the physics of superadiabatic ampli- 
fication ||. It is this phenomenon's vast enhancement of the usual 0-point 
energy which causes one of the zeroth order infrared solutions to dominate 
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at late times. These approximations therefore isolate precisely the leading 
late time infrared effect we wish to study. In fact this is all that can be 
reliably studied using quantum general relativity. The ultraviolet regime, 
which these approximations fail to capture, cannot in any case be described 
perturbatively by quantum general relativity. 

What remains is to implement the infrared approximations (|160| . |162| , |167p 



in the various propagator mode sums derived in Section 4. Since we are 
only computing the amputated 1-point function to one loop order, these 
propagators are all coincident. They may, however, bear derivatives. Since 
space derivatives add factors of k, which are small in the infrared, we need 
only consider time derivatives. We shall therefore set Ax = but keep 
the two times nonzero. With these conventions all the propagators can be 
described in the standard form: 

zA (x; x') — f(t) ■ g{t ! ) ■ J dk^ ■ h(k) , (168) 

where it should be noted that ii* and if* are functions of the co-moving wave 
number k, determined by the horizon crossing condition (|158|) . 

From ( |160[ ) and the mode sums ( |9"3"| , |8~iD which define it we see that the 



infrared part of iAa(x : x') approaches a constant: 



iA A (x; x') — > 1 ■ 1 ■ ^ J dk^j- ■ 1 . (169) 



The behavior of the B mode ( [L62|) and the mode sums (|9"3|J8~4] ) which define 
2 A# (a;; x') show that its infrared part actually falls off: 

iA B (x; x') — ► , } -. . • ; I - . ■ At / dk^l ■ k 2 . (170) 
v ' ' a {t)H(t) a Q (t')H(t') 2n 2 J 2k K ' 

Since the momentum integral is dominated by the ultraviolet, rather than 
the infrared, we conclude that the infrared part of this propagator is zero. 
The -000 propagator involves B modes (|162|) and C modes ( |167|) . From its 



defining relations (0, |4]) we determine its infrared part to be: 



• a / i\ -Hit) -Hit') 1 f 17 Hl Hi .^^ 
iAJx; x') — > — the ? / dk— 1- . (171) 
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The mixed propagator is defined by relations ( ]91"| , |8"4| ) . Again applying the 
infrared mode approximations ( |162| , |1671 ) gives: 



-H(t) J-H(V) 1 t Hi Hi 

lA ^ x) -> Jm • m~ -2^l dk 2k-^- (172) 

The same infrared limits, applied to its defining relations (|92"| , |34D , reduce the 
scalar propagator to: 



zA ( x - *') v ~ m . V- g ^ . J_ / dk S. . IL (173) 
ZA7(X ' Xj H{t) H(t') 2^J d 2k - H / { } 

By itself it is the strongest but one must allow for the effect of factors and 
derivatives from the interaction vertex. 

Since modes only become infrared after horizon crossing, the momentum 
integrations are cut off at k — H(t)a (t). The integral can be evaluated by 
first changing variables from k to the horizon crossing time £*: 

dk « H 2 {U)a (QdU , (174) 
and then employing the slow roll expansions of Section 2: 

4(4) « ±r*.*&, cm) 



2vr 2 7 2A; ' 4vr 2 Jo -H 

1 fHf-H 6 (t) 
24tt 2 \ ^2 



(176) 



Note that if is approximately constant for the quadratic potential considered 
by Mukhanov, Abramo and Brandenberger H. 



7 Amputated 1-point functions 

The purpose of this section is first to obtain one loop results for the three 



amputated 1-point functions (108- 109 ) defined in Section 5. We then ex- 



ploit the technology of Section 5 to compute the two observables HeS(t) 
( |106|) and $(t) ( |107| ). We begin by explaining how cubic interactions are used 
to compute the amputated 1-point functions. 
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At one loop order the amputated 1-point functions consist basically of 
coincident propagators contracted into cubic interaction vertices. In addition 
to the usual factor i there is an % from the kinetic operator acting on the 
external propagator. There is also an extra factor of k from the fact that we 
define the 1-point functions ( |102| - |103|) as (0|m/v|0) an d (O|k0|O). 



As an example let us consider the interaction — jKm 2 fl 4l (j) 2 ilj, which is one 
of the many cubic terms descending from the scalar mass. The external line 
can attach to any of the three quantum fields: c^^^rf 1 ' . When it attaches 
to one of the two scalar fields S(t) receives the following contribution: 

i ■ k ■ - l -Km 2 Vt 4 ■ 2 • (0|T [0(x)^(x')] |0) 

= ^K 2 m 2 nHAp(x; + ^)rr} ^ , (177) 

= {K 2 m 2 nHA,3(x;x')} x/ _ x . (178) 

Most of the coincindence limit is ultraviolet nonsense which quantum general 
relativity cannot be trusted to treat correctly and which must in any case 
have been subtracted off in order for inflation to begin in the first place. The 
time dependent, physically significant part comes from the superadiabatically 
amplified infrared modes. From equation ( |1 72] ) we see that the leading effect 
from these is: 

K W S A e( ,;,) -3*^(0 • tSJ • JL p^M) , 

k 2 [Hf-H 6 (t)\ ( 



?7T 2 V H 2 (t) 




4(t) . (179) 



In the notation used at the end of Section 5 there are contributions for 
N = -2 and N = +4 with coefficients of k 2 /(8tt 2 ) times +Hf and -1, 
respectively. 

When the external leg attaches to the pseudo-graviton field the interaction 
makes contributions to a{t) and 7(£). Because ip = — ip 00 + ipa these have 
opposite signs. The contribution for a{t) is: 

\i- k- --Km 2 Q 4 ■ (0\T[<f){x)<j){x')]\0)\ =-K 2 m 2 QHAJx;x), (180) 
I 4 J x ' =x 4 
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Table 1: Vertex factors contracted into ^axB^ai^oaBz with V'ai/?! external. 
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^-l^-m-M^ 1 )' (181) 

- «" ff-^Ut). (182) 



32tt 2 V # 2 (t) 

In the notation used at the end of Section 5 there are contributions for 
N = —2 and N = +4 with coefficients of k 2 /(32tt 2 ) times +Hf and -1, 
respectively. 

Tables 1-3 give the various cubic interaction vertices. Note that those of 
Table 1 have been partially symmetrized by making the external leg attach to 
pseudo-graviton #1. It should also be noted that derivatives with respect to 
the attached field are interpreted through integration by parts (IBP) as acting 
on the entire result. Further, spatial translation invariance (STI) allows the 
result to depend only upon the conformal time after the expectation value is 
taken. For example, in vertex # 3 of Table 1, the derivative with respect to 
line #1 is interpreted as follows: 

d a 2 ibp _d sri^d 
ox a2 dr\ 

Because of the r] ai131 vertex #3 contributes to a(t) and j(t) with opposite 
signs. The contribution to a(t) is: 

A {< • k ■ % -kH& • (o|t [Vw 2 WVwa^')] \o)t a n^ v a ^ 

= T~ {—i^HD, 3 ■ 2iA a (x; x')\ , (184) 
arj I 2 J x ' =x 

- «' W-^Ut). (186, 



8tt 2 V H 2 (t) 

When the contributions from all vertices are summed the amputated 1- 
point functions are: 

j* 3k 2 ( H*-H*(t) \ 4 
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Table 2: Vertex factors contracted into ip^^uo^uo, 
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Table 3: Cubic interactions involving 0. 
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32vr 2 



H 2 (t) 
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S(t) 



IR 



k 2 (Hf-H 6 (ty 



-H 



8vr 2 



H*{t) 



H(t) 



a$(t) . 



:i88) 



;i89) 



In the notation at the end of Section 5 this corresponds to the following 
coefficients: 

* 3k 2 



Ot-2 = -7-2 = ^5-2 



a+4 



-7+4 = ^5+4 



32vr 2 7 ' 
3k 2 



32tt 



2 ' 



(190) 
(191) 



From (|141[ - |142 ) we see that the cosmological expansion rate and the quantum- 
corrected scalar are: 



H cS (t) = H(t)ll 



72tt 2 



576vr 2 



H[ -3HfH 3 (t) + 2H 6 (ty 
f H[ -2HfH 3 (t) + H 6 (tY 



, (192) 
. (193) 



8 Summary and discussion 

This paper is first of all a check on the calculation of Mukhanov, Abramo 
and Brandenberger. We certainly agree with the sign and the leading time 
dependence that can be inferred for H e g(t) and from their published 
results 0, |3|. From unpublished work we see that even the numerical factors 
agree. It is worth emphasizing that we employed the standard formalism 
of covariant quantization while they used a truncated version of the canon- 
ical formalism. The two calculations were also done in completely different 
gauges: we added a covariant gauge fixing term whereas they used a physical 
gauge. It would be difficult to imagine two more completely different calcu- 
lational schemes. Yet we got the same results in the end, for both H e g(t) and 
$(t). This is an enormously powerful check on the validity of their 
work and on the physical reality of the effect. If this back-reaction is 
a gauge chimera it is a remarkably consistent one. 

Our motivation for this work was the excellent questions posed earlier this 
year by Unruh Although our calculation is itself a sort of answer we have 
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analyzed some of his arguments more generally in Section 3. In particular, it 
turns out that the seeming disagreement between Unruh's long wavelength 
solutions and those of Mukhanov et al. derives from different definitions for 
what is "zeroth order" when expanding in powers of the wave number. If one 
keeps higher terms in k 2 it turns out that Unruh's second solution implies 
an order k 2 result for the variable used by Mukhanov et al., and this result 
is just 4k 2 / k 2 times their second solution. So neither of the long wavelength 
solutions of Mukhanov et al. is unphysical. 

It is important to understand that a compelling physical mechanism un- 
derlies the back-reaction of Mukhanov et al.. It is the self-gravitation be- 
tween superadiabatically amplified long wavelength modes. A simple physi- 
cal model is that virtual particles whose wavelengths are comparable to the 
Hubble radius become trapped in the expansion of spacetime and are not 
able to recombine. As the particles are pulled apart their long range grav- 
itational potentials fill the intervening space, adding with the potentials of 
earlier pairs. Because gravitation is attractive these potentials resist the fur- 
ther expansion of spacetime, thereby slowing inflation. There is absolutely 
no question that this process should occur for quanta, such as gravitons 
and minimally coupled scalars, which lack conformal invariance but are still 
massless on the Hubble scale. The only issues concern the strength of the 
back-reaction, its time dependence, and whether or not it can eventually stop 
inflation. 

The analogous back-reaction has already been demonstrated for gravitons 
when inflation is driven by a positive bare cosmological constant M. In this 
case it comes at two loops because the pair creation event is already one loop 
and the absence of linearized mixing between the dynamical spin 2 gravitons 
and the spin gravitational potentials postpones self-gravitation to next 
order. When the superadiabatically amplified quanta are themselves scalar 
their mixing with the spin gravitational potentials allows self-gravitation 
to occur at one loop order. However, there does not have to be such a one 
loop effect ||18|| . The feature which seems to distinguish those scalar-driven 
models which show slowing at one loop from those which slow at two loops is 
the rate at which superadiabatic amplification injects 0-point energy. If this 
is less than or equal to the physical 3- volume's inflation then there is no one 
loop effect; if superadiabatic amplification injects 0-point energy faster than 
the 3- volume inflates to absorb it then there is a one loop effect. 

Of course gravitons presumably drive a two loop effect in this model as 
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well. There may also be significant scalar effects at higher loops. Higher loop 
processes are interesting in that they derive from the coherent superposition 
of interactions over the invariant volume of the past lightcone, which can 
grow arbitrarily large. There is no barrier to considering such questions in 
the covariant formalism we have developed. The formalism of Mukhanov et 
al. would have to be extended to make this possible. 

It may be of general interest that we were able to extend the Feynman 
rules of Iliopoulos et al. so that they apply to an arbitrary scalar potential. 
This was the work of Section 4. Of course we can only express the propagators 
as mode sums, where even the mode functions remain to be determined. 
But we have shown in Section 6 how these mode functions can be usefully 
expanded, both in the ultraviolet and in the infrared. And the calculation is 
an explicit example of how interesting effects can be obtained. It should now 
be possible to re-do the two loop computation of Tsamis and Woodard [|7| 
for an arbitrary background. This should completely determine the effective 
field equations needed to evolve past the end of inflation to arbitrarily late 
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